#Figure A1

data<-read.dta13("Remaining Seized Replication Data.dta")

plot<-na.omit(data[,c("agg_implem_wo", "lag_runner",
                      "sum_un", "lag_ltroop", "lag_lobs", 
                      "lag_lpolice", "sum_shame",
                      "sum_coord", "sum_resource",
                      "all_words", "conflict_type",
                      "dead_log", "war_dur_months",
                      "infant_wb", "merged_gdp_growth",
                      "lag_constraint", "donor_implem",
                      "io_or_civil_sigs", "PTS_S", "v2x_cspart",
                      "external_miss"
)])
names(plot)<-c("Aggregate Implementation Rate",
               "Total UNSC Resolutions", 
               "Log of Total UN Personnel",
               "Log of UN Troops",
               "Log of UN Observers",
               "Log of UN Police",
               "Naming and Shaming",
               "Coordinating", "Resources",
               "All Approaches", "Conflict Type",
               "Log of War Casualties", 
               "War Duration (Months)",
               "Infant Mortality", "GDP Growth",
               "Executive Constraints",
               "Donor Implementation",
               "Civil Society or IO Signatories",
               "Political Terror Scale (State Dept.)",
               "Civil Society Strength", "Third Party Support")

corrplot(as.matrix(cor(plot)), diag=F, type="upper", method="color", tl.col="black")

#Figure A1
pdf("FigureA1.pdf")
corrplot(as.matrix(cor(plot)), diag=F, type="upper", method="color", tl.col="black")
dev.off()

rm(fun,plot,data)
